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Abstract Persistent high-energy emission of magnetars 
is produced by a plasma corona around the neutron star, 
with total energy output of ~ 10 36 erg/s. The corona 
forms as a result of occasional starquakes that twist the 
external magnetic field of the star and induce electric 
currents in the closed magnetosphere. Once twisted, the 
magnetosphere cannot untwist immediately because of 
its self-induction. The self-induction electric field lifts 
particles from the stellar surface, accelerates them, and 
initiates avalanches of pair creation in the magnetosphere. 
The created plasma corona maintains the electric cur- 
rent demanded by curl B and regulates the self-induction 
e.m.f. by screening. This corona persists in dynamic equi- 
librium: it is continually lost to the stellar surface on 
the light-crossing time ~ s and replenished with 

new particles. In essence, the twisted magnetosphere acts 
as an accelerator that converts the toroidal field energy 
to particle kinetic energy. The voltage along the mag- 
netic field lines is maintained near threshold for igni- 
tion of pair production, in the regime of self-organized 
criticality. The voltage is found to be about ~ 1 GeV 
which is in agreement with the observed dissipation rate 
~ 10 36 erg/s. The coronal particles impact the solid 
crust, knock out protons, and regulate the column den- 
sity of the hydrostatic atmosphere of the star. The tran- 
sition layer between the atmosphere and the corona is 
the likely source of the observed 100 keV emission from 
magnetars. The corona also emits curvature radiation 
up to 10 14 Hz and can supply the observed IR-optical 
luminosity. 
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1 Introduction 

At least 10% of all neutron stars are born as magne- 
tars, with magnetic fields B > 10 14 G. Their activity is 
powered by the decay of the ultrastrong field and lasts 
about 10 4 years. They are observed at this active stage 
as either Soft Gamma Repeaters (SGRs) or Anomalous 
X-ray Pulsars (AXPs) [T]. 

Besides the sporadic X-ray outbursts, a second, per- 
sistent, form of activity has been discovered by studying 
the emission spectra of magnetars. Until recently, the 
spectrum was known to have a thermal component with 
temperature k^T ~ 0.5 keV, interpreted as blackbody 
emission from the star's surface. The soft X-ray spec- 
trum also showed a tail at 2 — 10 keV with photon index 
r = 2 — 4. This deviation from pure surface emission 
already suggested that energy is partially released above 
the star's surface. Recently, observations by RXTE and 
INTEGRAL have revealed even more intriguing feature: 
magnetars are bright, persistent sources of 100 keV X- 
rays PJ. This high-energy emission forms a separate com- 
ponent of the magnetar spectrum. It becomes dominant 
above 10 keV, has a very hard spectrum, r ~ 1, and 
peaks above 100 keV where the spectrum is unknown. Its 
luminosity, L ~ 10 36 erg s , even exceeds the thermal 
luminosity from the star's surface. The observed hard X- 
rays can be emitted only in the exterior of the star and 
demonstrate the presence of an active plasma corona. 

This corona can form as a result of starquakes that 
shear the magnetar crust and its external magnetic field. 
The persistence of magnetospheric twists can explain 
several observed properties of magnetars. Thompson, Lyu- 
tikov, & Kulkarni investigated the observational con- 
sequences using a static, force-free model, idealizing the 
magnetosphere as a globally twisted dipole. They showed 
that a twist affects the spindown rate of the neutron star: 
it causes the magnetosphere to flare out slightly from a 
pure dipole configuration, thereby increasing the braking 
torque acting on the star. Although the calculations of 
force-free configurations are independent of the plasma 
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behavior in the magnetosphere, they rely on the pres- 
ence of plasma that can conduct the required current 
j B = (c/47r)V x B. This requires a minimum particle 
density n c = js/cc. Thompson et al. [3] showed that 
even this minimum density can modify the stellar radia- 
tion by multiple resonant scattering. 

The problem of plasma dynamics in the closed twisted 
magnetosphere has been formulated in jjj, and a simple 
solution has been found to this problem. It allows one 
to understand the observed energy output of the mag- 
netar corona and the evolution of magnetic twists. It is 
worth mentioning that plasma behavior around neutron 
stars has been studied for decades in the context of radio 
pulsars, and that problem remains unsolved. The princi- 
ple difference with radio pulsars is that their activity is 
caused by rotation, and dissipation occurs on open mag- 
netic lines that connect the star to its light cylinder. 
Electron-positron pairs are then created on open field 
lines By contrast, the formation of a corona around 
a magnetar does not depend on its rotation. All observed 
magnetars are slow rotators (J? = 2-k/P ~ 1 Hz), and 
their rotational energy is unable to feed the observed 
coronal emission. Practically the entire plasma corona 
is immersed in the closed magnetosphere and its heat- 
ing must be caused by some form of dissipation on the 
closed field lines. It is this closure that facilitates the 
corona problem by providing both boundary conditions 
at the two footpoints of a field line. 

The basic questions that one would like to answer are 
as follows. How is the magnetosphere populated with 
plasma? The neutron star surface has a temperature 
k-g,T Ss 1 keV and the scale-height of its atmospheric 
layer (if it exists) is only a few cm. How is the plasma sup- 
plied above this layer and what type of particles populate 
the corona? How is the corona heated, and what are the 
typical energies of the particles? If the corona conducts 
a current associated with a magnetic twist V x B / 0, 
how rapid is the dissipation of this current, i.e. what is 
the lifetime of the twist? Does its decay imply the dis- 
appearance of the corona? 

An outline of the proposed model is as follows. The 
key agent in corona formation is an electric field E\\ 
parallel to the magnetic field. En is generated by the 
self-induction of the gradually decaying current and in 
essence measures the rate of the decay. It determines the 
heating rate of the corona via Joule dissipation. If E\\ = 
then the corona is not heated and, being in contact with 
the cool stellar surface, it will have to condense to a thin 
surface layer with k^T 5s 1 keV. The current-carrying 
particles cannot flow upward against gravity unless a 
force eEu drives it. On the other hand, when En exceeds 
a critical value, avalanches are triggered in the mag- 
netosphere, and the created pairs screen the electric field. 
This leads to a "bottleneck" for the decay of a magnetic 
twist, which implies a slow decay. 

Maintenance of the corona and the slow decay of the 
magnetic twist are intimately related because both are 



governed by Eu . In order to find E\\ , one can use Gauss' 
law V • E = 4irp where p is the net charge density of the 
coronal plasma. This constraint implies that E and p 
must be found self-consistently. The problem turns out to 
be similar to the classical double-layer problem of plasma 
physics with a new essential ingredient: creation. 

A direct numerical experiment can be designed that 
simulates the creation and behavior of the plasma in the 
magnetosphere. The experiment shows how the plasma 
and electric field self-organize to maintain the time-average 
magnetospheric current j = j B demanded by the mag- 
netic field, (4%/c) j = V x B. The electric current admits 
no steady state on short timescales and keeps fluctuating, 
producing e^ avalanches. This state may be described as 
a self-organized criticality. Pair creation is found to pro- 
vide a robust mechanism for limiting the voltage along 
the magnetic lines to e<P e 1 GeV and regulate the 
observed luminosity of the corona. 



2 Mechanism of corona formation 

A tightly wound-up magnetic field is assumed to exist 
inside magnetars 0. The internal toroidal field can be 
much stronger than the external large-scale dipole com- 
ponent. The essence of magnetar activity is the trans- 
fer of magnetic helicity from the interior of the star to 
its exterior, where it dissipates. This involves rotational 
motions of the crust, which inevitably twist the exter- 
nal magnetosphere anchored to the stellar surface. The 
magnetosphere is probably twisted in a catastrophic way, 
when the internal magnetic field breaks the crust and 
rotates it. Such starquakes are associated with observed 
X-ray bursts [Hj- The most interesting effect of a star- 
quake for us here is the partial release of the winding 
of the internal magnetic field into the exterior, i.e., the 
injection of magnetic helicity into the magnetosphere. 

Since the magnetic field is energetically dominant in 
the magnetosphere, it must relax there to a force-free 
configuration with j x B = 0. Electric current j initi- 
ated by a starquake flows along the magnetic field lines 
to the exterior of the star, reaches the top of the field line, 
and comes back to the star at the other footpoint, then 
enters the atmosphere and the crust, following the mag- 
netic field lines ( j || B). Sufficiently deep in the crust the 
current flows across the field lines and changes direction 
so that it can flow back to the twisted footpoint through 
the magnetosphere (or directly through the star, depend- 
ing on the geometry of the twist). The current can flow 
across the field lines in the deep crust because it is able 
to sustain a significant Ampere force j x B/c ^ — this 
force is balanced by elastic forces of the deformed crust. 

The emerging currents are easily maintained during 
the X-ray outburst accompanying a starquake. A dense, 
thermalizcd plasma is then present in the magnetosphere, 
which easily conducts the current. Plasma remains sus- 
pended for some time after the starquake because of the 
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transient thermal afterglow 9 j . Eventually the afterglow 
extinguishes and the radiative flux becomes unable to 
support the plasma outside the star. The decreasing den- 
sity then threats the capability of the magnetosphere to 
conduct the current of the magnetic twist. A minimal 
"corotation" charge density p co = — ft ■ B/2nc is al- 
ways maintained |10| . but it is far smaller than needed 
to supply the current j ~ (c/4tt)(B/R^s)- Can the lack 
of charges stop the flow of electric current? A simple 
estimate shows that the current cannot stop quickly be- 
cause of its self-induction. A slow decay of the current 
generates a sufficient self-induction voltage that helps the 
magnetosphere to re-generate the plasma that carries the 
current, so that the twisted force- free configuration per- 
sists with a qausi-steady j ^ 0. 

The stored energy of non-potential (toroidal) mag- 
netic field associated with the ejected current is sub- 
ject to gradual dissipation. In our model, this dissipation 
feeds the observed activity of the corona. The stored en- 
ergy is concentrated near the star and carried by closed 
magnetic lines with a maximum extension radius i? max ~ 
2Rns- Thus, most of the twist energy will be released if 
these lines untwist, and we focus here on the near mag- 
netosphere i? max ~ 2Rns- 

Consider a magnetic flux tube with cross section S ^ 
R^ s and length L ~ _Rns which carries a current / = 
Sj. The stored magnetic energy of the current per unit 
length of the tube is 



f T 2 



s. 



(1) 



L c 2 

The decay of this energy is associated with an electric 
field parallel to the magnetic lines E||: this field can ac- 
celerate particles and convert the magnetic energy into 
plasma energy. Conservation of energy can be expressed 
as 

~j-V-(c^], (2) 



dt V 8tt 



E 



4tt 



as follows from Maxwell's equations with E <§C B. The 
first term on the right-hand side is the Joule dissipation 
caused by E|| . The second term — the divergence of the 
Poynting flux — vanishes for the induced field E||. 

The decay of the twist is related to the voltage be- 
tween the footpoints of a magnetic line, 

*e = / E dl. (3) 

J A 

Here A and C are the anode and cathode footpoints and 
dl is the line element; the integral is taken along the 
magnetic line outside the star (Fig. 1). The product <t> e I 
approximately represents the dissipation rate £ in the 
tube. (The two quantities are not exactly equal in a time- 
dependent circuit.) The instantaneous dissipation rate is 
given by 

£= / C E-I(0dl, (4) 




Fig. 1 Schematic picture of a current-carrying magnetic 
tube anchored to the star's surface. The current is initiated 
by a starquake that twists one (or both) footpoints of the 
tube. The current flows along the tube up to the magneto- 
sphere and comes back to the star at the other footpoint. 
A self-induction voltage is created along the tube between 
its footpoints, which accompanies the gradual decay of the 
current. The voltage is generated because the current has a 
tendency to become charge-starved above the atmospheric 
layer whose scale-height h — ksT ' j gm p is a few cm. 



where 1(1) — Sj is the instantaneous current at position 
I in the tube. The current is fixed at footpoints A and 
C, I = Jo; it may, however, fluctuate between the foot- 
points. 

The voltage <£ e is entirely maintained by the self- 
induction that accompanies the gradual untwisting of the 
magnetic field. Thus, <P e reflects the gradual decrease of 
the magnetic helicity in the tube, 
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/ E • B S(l) dl. 

J A 



(5) 



(Here dV = S(l)dl is the volume element.) It also de- 
termines the effective resistivity of the tube: 1Z = ^ e /I- 
The self-induction voltage passes the released magnetic 
energy to charged particles, and a higher rate of untwist- 
ing implies a higher energy e<£ e gained per particle. A 
huge magnetic energy is stored in the twisted tube, and 
a quick untwisting would lead to extremely high Lorentz 
factors of the accelerated particles. 

There is, however, a bottleneck that prevents a fast 
decay of the twist: the tube responds to high voltages 
through the copious production of pairs. Runaway 
pair creation is ignited when electrons are accelerated to 
a certain Lorentz factor j± ~ f0 3 (see below). The cre- 
ated e ± plasma does not tolerate large Eu — the plasma 
is immediately polarized, screening the electric field. This 
provides a negative feedback that limits the magnitude 
of <P e and buffers the decay of the twist. 

A minimum <P e is needed for the formation of a corona. 
Two mechanisms can supply plasma: (1) Ions and elec- 
trons can be lifted into the magnetosphere from the stel- 
lar surface. This requires a minimum voltage, 



e# e ~ grriiRns ~ 200 MeV, 
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corresponding to En that is strong enough to lift ions (of 
mass mi) from the anode footpoint and electrons from 
the cathode footpoint. (2) Pairs can be created in the 
magnetosphere if 



e<P e — 7 ± toc 2 = 0.57± MeV, 



(7) 



which can accelerate particles to the Lorentz factor y± 
sufficient to ignite e creation. 

If <P e is too low and the plasma is not supplied, a 
flux tube is guaranteed to generate a stronger electric 
field. The current becomes slightly charge-starved: that 
is, the net density of free charges becomes smaller than 
| V x B|/47re. The ultrastrong magnetic field, whose twist 
carries an enormous energy compared with the energy of 
the plasma, does not change and responds to the decrease 
of j by generating an electric field: a small reduction of 
the conduction current j induces a displacement current 
(l/47r)9E/<9i = (c/4?r)V x B — j. The longitudinal elec- 
tric field Eu then quickly grows until it can pull particles 
away from the stellar surface and ignite pair creation. 

The limiting cases En — > (no decay of the twist) 
and -Eii — » oo (fast decay) both imply a contradiction. 
The electric field and the plasma content of the corona 
must regulate each other toward a self-consistent state, 
and the gradual decay of the twist proceeds through a 
delicate balance: E\\ must be strong enough to supply 
plasma and maintain the current in the corona; how- 
ever, if plasma is oversupplied En will be reduced by 
screening. A cyclic behavior is possible, in which plasma 
is periodically oversupplied and En is screened. Our nu- 
merical experiment shows that such a cyclic behavior 
indeed takes place. 



3 Electric circuit: numerical experiment 

Three facts facilitate the simulation of the coronal elec- 
tric circuit: 

1. The ultrastrong magnetic field makes the parti- 
cle dynamics 1-D. The magnetic field lines are not per- 
turbed significantly by the plasma inertia, and they can 
be thought of as fixed "rails" along which the particles 
move. The particle motion is confined to the lowest Lan- 
dau level and is strictly parallel to the field (the lifetime 
of a particle in an excited Landau state is tiny). 

2. The particle motion is collisionless in the magneto- 
sphere. It is governed by two forces only: the component 
of gravity projected onto the magnetic field and a collec- 
tive electric field E\\ which is determined by the charge 
density distribution and must be found self-consistently. 

3. The star possesses a dense and thin atmospheric 
layer. 1 Near the base of the atmosphere, the required cur- 
rent j B = (c/47r)V x B is easily conducted, with almost 

1 Such a layer initially exists on the surface of a young 
neutron star. Its maintenance is discussed in 



no electric field. Therefore the circuit has simple bound- 
ary conditions: E\\ — and fixed current. The atmo- 
sphere is much thicker than the skin depth of the plasma 
and screens the magnetospheric electric field from the 
star. 

Our goal is to understand the plasma behavior above 
the screening layer, where the atmospheric density is ex- 
ponentially reduced and an electric field En must de- 
velop. The induced electric field E and conduction cur- 
rent j satisfy the Maxwell equation, 

^ 4tt . 1 <9E 
VxB-— j + - — . 8 

c c at 

Here j is parallel to the direction of the magnetic field, 
and the force- free condition requires V x B || B. There- 
fore, in the fixed magnetic configuration where V x B 
does not vary with time, d E/dt is parallel to B, i.e. only 
Eu is created by the self-iduction effect. We consider a 
simple model problem where no pre-existing perpendicu- 
lar field E^ exists; more precisely, we require V- Ex — 0. 
Then Gauss' law reads 

47rp = V- E|, = -JL, (9) 



dl 1 

where I is length measured along the magnetic tube. 
Then the problem becomes strictly 1-D since Ej_ has no 
relation to charge density and falls out from the prob- 
lem. Note that the approximation of ID circuit, jsit) = 
const, where jb = (c/47r)|V x B|, excludes the excita- 
tion of transverse waves in the magnetosphere, which in 
reality can exist. These waves are described by the cou- 
pled fluctuations 8js and Ej_ on scales much smaller 
than the circuit size L. (Wavelengths A ~ c/ojp can be 
excited by plasma oscillations.) These fluctuations may 
be expected to have a small effect on the circuit solution 
unless 5jb becomes comparable to jg. 

With js(t) — const, the particle motions on different 
magnetic lines are decoupled. Indeed, the particle dy- 
namics is controlled by electric field E = Eu which is re- 
lated to the instantaneous charge distribution by Gauss' 
law ©. Thus, E(l) on a given magnetic line is fully de- 
termined by charge density p{l) on the same line. The 
plasma and electric field evolve along the line as if the 
world were 1-D. 

If the conduction current j is smaller than js, then 
dE\\/dt > and an electric field appears that tends to 
increase the current. Alternatively, if j > js then an 
electric field of the opposite sign develops which tends 
to reduce the current. Thus, j is always regulated to- 
ward 2b- This is the standard self-induction effect. The 
timescale of the regulation, r, is very short, of the order 



of the Langmuir oscillation timescale ui } 
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cal deviations from charge neutrality tend to be erased 
on the same plasma timescale. The net charge imbalance 
of the corona must be extremely small — a tiny imbal- 
ance would create an electric field that easily pulls out 
the missing charges from the surface. When the corota- 
tion charge density is neglected, the net charge of the 
corona vanishes. 
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Fig. 2 Set up of numerical experiment. Thin and dense 
plasma layers are maintained near the cathode and anode by 
injecting cold particles through the boundaries of the tube 
(footpoints of a magnetic flux tube). The electric current is 
kept constant at the boundaries and the system is allowed 
to evolve in time until a quasi-steady state is reached. With- 
out voltage between anode and cathode, the current cannot 
flow because gravity g traps the particles near the bound- 
aries. The constant current at the boundaries, however, im- 
plies that voltage is immediately generated should the flow 
of charge stop in the tube. The experiment aims to find the 
self-induction voltage that keeps the current flowing. 



Adopting a characteristic velocity v ~ c of the coro- 
nal particles, the Debye length of the plasma Ad — v/lup 
can be taken equal to the plasma skin depth c/uop. An 
important physical parameter of the corona is the ratio 
of the Debye length to the circuit size L ~ i?NS ; 



c = 



L 



Luip 



< 1. 



(10) 



The force-free condition j B x B together with V • 
j B = requires that jp(l) oc B(l) along the magnetic 
line. We can scale out this variation simply dividing all 
local quantities (charge density, current density, and elec- 
tric field) by B. This reduces the problem to an equiv- 
alent problem where js(0 = const. Furthermore, only 
forces along magnetic lines control the plasma dynam- 
ics, and the curvature of magnetic lines falls out from 
the problem. Therefore, we can set up the experiment so 
that plasma particles move along a straight line connect- 
ing anode and cathode (Fig. 2). We designated this line 
as the z-axis, so that I = z. 

This one-dimensional system may be simulated nu- 
merically. Consider N particles moving along the z-axis 
(N ~ 10 6 in our simulations). The electric field acting 
on a particle at a given z is given by 



E{z) = 4tt / 
Jo 



pdz, 



(11) 



where E(0) = E{L) = are the boundary conditions. 
Given the instantaneous positions Zi of all charges ei we 
immediately find E(z). To find a quasi-steady state of 




Fig. 3 Circuit without e production. Upper panel shows 
the normalized charge density and lower panel shows the elec- 
tric field in units of E g defined by eE g = rriigo. Solid and 
dotted curves correspond to two different moments of time, 
demonstrating that a quasi-steady state has been reached. In 
this simulation, rrii = 10m e and C, = c/ujpL ~ 0.01. The final 
state is the relativistic double layer described analytically in 

rm 



such a system we let it relax by following the particle 
motion in the self-consistent electric field. 

Consider first a circuit where e^- production is (ar- 
tificially) forbidden, so that only the cold hydrostatic 
atmosphere can supply particles to the corona. We find 
that such a circuit quickly relaxes to a state where it 
acts as an ultra-relativistic linear accelerator (Fig. 3). 
In this state, oscillations of electric field (on the plasma 
timescale uip 1 ) are confined to the thin atmospheric lay- 
ers. A static accelerating electric field is created above 
the layers where the atmosphere density is exponentially 
reduced. The asymmetry of the solution is caused by the 
difference between the electron and ion masses. 

This configuration is a relativistic double layer [IT] . 
It is well described by Carlqvist's solution, which has no 
gravity in the circuit and assumes zero temperature at 
the boundaries, so that particles are injected with zero 
velocity. According to this solution, the potential drop 
between anode and cathode, in the limit e<P e 3> m^c 2 , is 



Zm, 



1/2 



1 



Lloi 



■ m e c, uj p 



(12) 



where Z is the ion charge number (in our experiment 
Z = 1 was assumed). The established voltage is much 
larger than is needed to overcome the gravitational bar- 
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rier <P g . It does not even depend on <P g as long as <P g is 
large enough to prohibit the thermally-fed regime. Grav- 
ity causes the transition to the linear-accelerator state, 
but the state itself does not depend on <P g . 

Thus, the first result of numerical experiment is that 
an ion-electron circuit (without e ± creation) in a gravi- 
tational field relaxes to the double layer of macroscopic 
size L and huge voltage <P e . The system does not find 
any state with a lower <£ e > even though it is allowed to 
be time-dependent and strong fluctuations persist in the 
atmospheric layers (reversing the sign of E). 

If the linear accelerator were maintained in the twisted 
magnetosphere, the twist would be immediately killed 
off: the huge voltage implies a large untwisting rate (§ 2). 
The electron Lorentz factor developed in the linear ac- 
celerator is (taking the real rrii/m e — 1836 and ( ~ 
3 x 1(T 9 ), 



7e 



e<Z> e 



20 



6 x 10 J 



(13) 



However, new processes will become important before 
the particles could acquire such high energies: produc- 
tion of pairs will take place. Therefore, the linear- 
accelerator solution cannot describe a real magnetosphere. 
We conclude that pair creation is a key ingredient of the 
circuit that will regulate the voltage to a smaller value. 



4 Pair discharge and self-organized criticality 




time 




Fig. 4 Spacetime diagrams illustrating the critical character 
of the e breakdown and formation of avalanches. The seed 
electron is placed at cathode and is accelerated to 7 ros after 



passing the distance a T , 



(7rcs — l)m e c jeE. The figure 



shows the worldlines of created particles, (a) a TCS /L = 0.35 
(supercritical case), (b) a les /L = 0.6 (critical case). 



In the magnetospheres of canonical radio pulsars with 
B ~ 10 12 G, e ± pairs are created when seed electrons 
are accelerated to large Lorentz factors 7 e ~ 10 7 . Such 
electrons emit curvature gamma rays that can convert 
to e ± off the magnetic field. In stronger fields, another 
channel of creation appears. It is also two-step: an 
accelerated particle resonantly upscatters a thermal X- 
ray photon, which subsequently converts to a pair 
0] . The resonant scattering requires the electron to have 
a Lorentz factor, 



7re 



B/B QED m e c 2 
1 - cos 6 kB huj x 



10 3 B u 



10 keV 



x 



(14) 



where Tilox and 6kB are the energy and pitch angle of the 
target photon. The distance over which an accelerated 
electron creates a pair is approximately equal to the free 
path to scattering l rcs . It depends only on B and the 
spectrum of target photons; Z res <C -Rns for the relevant 
parameter Bi 5 hivx < 100 keV Q]. 

When the voltage in the circuit becomes high enough 
to accelerate electrons to 7 res , an e ± breakdown develops 
in the magnetic flux tube. An illustrative toy model of 
breakdown is shown in a spacetime diagram in Figure 4. 
Each pair-creation event gives two new particles of op- 
posite charge, which initially move in the same direction. 
One of them is accelerated by the electric field and lost 
after a time < L/c, while the other is decelerated and can 



reverse direction before reaching the boundary. This re- 
versal of particles in the tube and repeated pair creation 
allows the e ± plasma to be continually replenished. In 
the super-critical regime (left panel in Fig. 4) more than 
one reversing particle is created per passage time L/c 
and an avalanche develops exponentially on a timescale 
~ a ros /c. In the near-critical regime shown in the right 
panel, just one reversing particle is created per passage 
time L/c. This critical state is unstable: sooner or later 
the avalanche will be extinguished. 

This toy model shows an essential property of the 
breakdown: it is a critical stochastic phenomenon. 
Above a critical voltage pair creation proceeds in a run- 
way manner, and the current and the dissipation rate 
would run away if the voltage were fixed. Below the crit- 
ical voltage pair creation does not ignite. The criticality 
parameter is L/a ICS — e$ e / (ires — l)m e c 2 . The tube with 
enforced current at the boundaries must self-organize to 
create pairs in the near-critical regime e<P e ~ 7 res m e c 2 
and maintain the current. 

The discharging tube is similar to other phenomena 
that show self-organized criticality, e.g., a pile of sand on 
a table If sand is steadily added, a quasi-steady state 
is established with a characteristic mean slope of the 
pile. The sand is lost (falls from the table) intermittently, 
through avalanches — a sort of "sand discharge." In our 
case, charges of the opposite signs are added steadily 
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instead of sand (fixed j at the boundaries), and voltage 
<P e = EL plays the role of the mean slope of a pile. 
The behavior of the discharging system is expected to 
be time-dependent, with stochastic avalanches. 

We implemented the process of pair creation in our 
numerical experiment. Regardless the details of this pro- 
cess (values of 7 res and l Ies ) the circuit relaxed to the 
state of self-organized criticality after a few t<j yn = L jc. 
This state is time-dependent on short timescales, but it 
has a well-defined steady voltage when averaged over a 
few t dyn (Fig. 5), 

e^ e ~ 7 rcs TO e c 2 . (15) 

We found that this relation applies even to circuits with 
7 res m e c 2 3> rriiC 2 , where lifting of the ions is energet- 
ically preferable to pair creation. We conclude that the 
voltage along the magnetic tube is self-regulated to a value 
just enough to maintain pair production and feed the cur- 
rent with e pairs. In all cases, the pair creation rate 
2N + ~ 3b I & is maintained. 

The solution for the plasma dynamics in the corona is 
strongly non-linear and time-dependent. It is essentially 
global in the sense that the plasma behavior near one 
footpoint of a magnetospheric field line is coupled to the 
behavior near the other footpoint. Remarkably, this com- 
plicated global behavior may be described as essentially 
one-dimensional electric circuit that is subject to simple 
boundary conditions E « on the surface of the star 
and can be studied using a direct numerical experiment. 

The current is carried largely by e ± everywhere in 
the tube. The ion fraction in the current depends on 
the ratio e<P e /miC 2 ~ r ) Tes m e /mi (Fig. 6). If this ratio 
is small, pairs are easily produced with a small electric 
field, too small to lift the ions, and the ion current is 
suppressed. In the opposite case, ions carry about 1/2 of 
the current. By coincidence, 7 ros mi/m e is typical for 
magnetars, which implies that ions carry ^ 10% of the 
current. 



5 Implications 

The basic finding of this work is that an e^ 1 corona must 
be maintained around a magnetar. It exists in a state 
of self-organized criticality, near the threshold for 
breakdown. The stochastic e discharge continually re- 
plenishes the coronal plasma, which is lost on the dy- 
namic time tdyn ~ 10~ 4 s. The established voltage along 
a twisted magnetic tube is marginally sufficient to ac- 
celerate an electron (or positron) to the energy where 
collisions with ambient X-ray photons spawn new pairs. 
This physical condition determines the rate of energy dis- 
sipation in the corona. The generation of this voltage is 
precisely the self-induction effect of the gradually decay- 
ing magnetic twist, and the energy release in the corona 
is fed by the magnetic energy of the twisted field. 

The rate of energy dissipation in the twisted mag- 
netosphere is given by Ldiss = I&e where I is the net 




time [L/ c] 



Fig. 5 Time history of one numerical experiment. The cir- 
cuit parameters are: mi — 10m e , 7 ros = mi/m e , £ = 0.01. 
After a few dynamical times the voltage <& e stops growing 
and the circuit enters an oscillatory regime. During each os- 
cillation, an increased voltage triggers e^ 1 discharge, then N+ 
drops, and <P e begins to grow again. The quasi-steady os- 
cillating state (self-organized criticality) is established. The 
horizontal lines show the voltage e^ e = 7rcs"ieC 2 and the 
minimum pair production rate 2N+ = jg /e that can feed the 
required current js- 

current through the corona. The current may be esti- 
mated as I ~ jbo? ~ Acj) (B/Rns) a 2 , where a is 
the size of a twisted region on the stellar surface and 
A<p = |V x B|(S/i?Ns) _1 ~ 1 characterizes the strength 
of the twist. The calculated e^ e ~ 1 GeV implies 

(*> 

Observed luminosity Ldi ss ~ 10 36 erg s _1 is consistent 
with a partially twisted magnetosphere, a ~ 0.3i?NSj or 
a global moderate twist with Acj) (e<P e /GeV) ~ 0.1. 

Once created, a magnetospheric twist has a relatively 
long but finite lifetime. The energy stored in it is £twist ~ 
(7 2 /c 2 )i?NS- This energy dissipates in a time 

Way = Id" * 3 { 10- ergs- J loevj y " 

Note that a stronger twist (brighter corona) lives longer. 
If the time between large-scale starquakes is longer than 
£decay> the magnetar should be seen to enter a quiescent 
state. Such behavior has been observed in AXP J1810- 
197: an outburst was followed by the gradual decay on a 
year timescale [Tl] . 
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Fig. 6 Fraction of the electric current carried by ions vs. 
voltage. In the left panel voltage is taken in units of rriiC 2 , 
and in the right panel — in units of 7 rcs m e c 2 . Different colors 
and symbols show circuits with nrii/m e — 10,30 and various 
values of 7 rcs ranging from m,/5m e to 5mi/m e (see [1]). 



The voltage <P e implies a certain effective resistivity of 
the corona, 1Z = <P e /I- This resistivity leads to spreading 
of the electric current across the magnetic lines, which is 
described by the induction equation dH/dt = — cV x E. 
The timescale of twist spreading to the light cylinder is 
comparable to idccay Therefore, the impact of a star- 
quake on spindown may appear with a delay of ~ years. 

Even a small ion current implies a large mass trans- 
fer through the magnetosphere over the ~ 10 4 -yr active 
lifetime of a magnetar. It requires a significant excava- 
tion of the crust at the anode footpoints of twisted mag- 
netic tubes. The bombarding relativistic electrons from 
the corona cause this excavation. They spall heavy ions 
in the uppermost crust, re-generate the light-element at- 
mosphere on the surface, and regulate its column density 
to a value ~ 100m p /aT 09- 

The radiative output from the corona likely peaks in 
its inner region because the coronal current is concen- 
trated on closed field lines with a maximum extension 



< 



2-Rns (§ 2). Emission from the inner corona must 



be suppressed above ~ 1 MeV, regardless the mecha- 
nism of emission, because photons with energy <; 1 MeV 
cannot escape the ultra-strong magnetic field. The ob- 
served high-energy radiation extends above 100 keV 
I15U16II17II18| . There is an indication for a cutoff between 
200 keV and 1 MeV from COMPTEL upper limits for 
AXP 4U 0142+61. 

The possible mechanisms of emission are strongly 
constrained. The corona has a low density n ~ n c = 



collisions are rare, so two-body radiative processes are 
negligible. The upscattering of keV photons streaming 
from the surface does not explain the observed 100-keV 
emission 0]. A possible source of the 100 keV X-rays 
is the transition layer between the corona and the ther- 
mal photosphere [T§1I3| . It is heated by the coronal beam 
and cooled by heat conduction toward the surface, and 
its temperature is regulated to ~ 100(£/£coui)~ 2 ^ 5 keV, 
where t/lcou\ < 1 parameterizes the suppression of ther- 
mal conductivity by plasma turbulence. The emitted X- 
ray spectrum may be approximated as single-temperature 
optically thin bremsstrahlung Its photon index be- 
low the exponential cutoff is close to —1, in agreement 
with spectra observed with INTEGRAL and RXTE. 

The observed pulsed fraction increases toward the 
high-energy end of the spectrum and approaches 100% at 
100 keV (Kuiper et al. 2004). If the hard X-rays are pro- 
duced by one or two twisted spots on the star surface, the 
large pulsed fraction implies that the spots almost dis- 
appear during some phase of rotation. Then they should 
not be too far from each other. They cannot be, e.g., 
antipodal because a large fraction of the star surface is 
visible to observer due to the gravitational bending of 
light: Sria/471-iZgra = (2 - 4GA//c 2 i?Ns) _1 « 3/4 [21]. 

The transition layer at a footprint of a coronal flux 
tube may reach high enough temperature for thermal 

production @|. The created pairs then elevate above 
the surface and can conduct the magnetospheric current. 
This opens a possibility for the corona to be fed by sur- 
face pair creation instead of the discharge. 

The observed infrared (i\-band) and optical lumi- 
nosities of magnetars are ~ 10 32 erg/s |21) . which is 
far above the Rayleigh-Jeans tail of the surface black- 
body radiation. The inferred brightness temperature is ~ 
10 13 K if the radiation is emitted from the surface of the 
star. Four principle emission mechanisms may then be 
considered: coherent plasma emission, synchrotron emis- 
sion by electrons, cyclotron emission by ions, and curva- 
ture radiation by bunches in the corona. The last two 
are most plausible (see @] for detailed discussion) . In par- 
ticular, curvature radiation in the inner magnetosphere 
may extend into the optical-IR range. Charges moving 
on field lines with a curvature radius Rc Ss i?NS will emit 
radiation with frequency vc ~ j^c/2ttRc, which is in the 
K-band [y = 10 14 Hz) or optical when 7 e m e c 2 ij GeV. 
The radio luminosity of a normal pulsar can approach 
~ 10~ 2 of its spindown power, and it is not implausible 
that the observed fraction ~ 10~ 4 of the power dissipated 
in a magnetar corona is radiated in optical-IR photons 
by the same coherent mechanism. 

The corona luminosity is simply proportional to the 
net current flowing through the magnetosphere, and is 
determined by the competition between occasional twist- 
ing (due to starquakes) and gradual dissipative untwist- 
ing of the magnetosphere. During periods of high ac- 
tivity, when starquakes occur frequently, the magneto- 
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spheric twist may grow to the point of a global insta- 
bility: when a critical A<p ?Z 1 is achieved, the magne- 
tosphcre suddenly relaxes to a smaller twist angle. A 
huge release of energy must then occur, producing a gi- 
ant flare; a model of how this can happen is discussed in 

m 

A possible way to probe the twist evolution is to mea- 
sure the history of spindown rate P of the magnetar. The 
existing data confirm the theoretical expectations: the 
spin-down of SGRs was observed to accelerate months to 
years following periods of activity |23U24j . A similar effect 
was also observed following flux enhancement in AXPs 
|25|. This "hysteresis" behavior of the spindown rate may 
represent the delay with which the twist is spreading to 
the outer magnetosphere. 

Since the density of the plasma corona is proportional 
to the current that flows through it, changes in the twist 
also lead to changes in the coronal opacity. The density 
of the plasma corona is close to its minimum n c = j / 'ec, 
and it is largely made of pairs rather than electron-ion 
plasma. The opacity should increase following bursts of 
activity and affect the X-ray pulse profile and spectrum. 
Such changes, which persist for months to years, have 
been observed following X-ray outbursts from two SGRs 

1251 ■ 
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